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Abstract 

We study the finite-size effects in some scaling systems, and sliow that the finite 
number of agents N leads to a cut-off in the upper value of the Pareto law for the 
relative individual wealth. The exponent a of the Pareto law obtained in stochastic 
multiplicative market models is crucially affected by the fact that N is always finite 
in real systems. We show that any finite value of A'^ leads to properties which can 
differ crucially from the naive theoretical results obtained by assuming an infinite 
N. In particular, finite N may cause in the absence of an appropriate social policy 
extreme wealth inequality a < 1 and market instability. 
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Power laws appear ubiquitously as probability distributions in a wide range of 
systems. They have been studied for a long time both theoretically and prac- 
tically. For instance, the relative number of people having a wealth between w 
and w + dw has been found already 100 years ago to fulfill the Pareto power 
law [1]: 

P{w)dw ~ w-^-''dw. (1) 



From the very beginning it turned out that especially for small values of 
exponents {a < 2) which imply slow decay of the distribution at infinity, this 
distribution form presents quite a few surprises, paradoxes and unexpected 
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properties. In particular, the limit of an infinite number of individuals — > oo 
is not defined for many of the measurements. 

In the present paper we study some dynamical processes where the naive 
theoretical results obtained for = oo fail for any arbitrarily large finite A^. 
One of the main sources of this non-uniform limit is the upper cut-off effect 
which any finite value of implies. We mean by "the upper cut-off effect" 
the following simple fact: Suppose the wealth of individuals in a particular 
system is distributed by a power law with o; < 2. Then, for any value w, there 
is a finite probability that an individual i possesses a wealth Wi larger than 
w. For a finite system, if one defines by Xi — Wi/W the relative wealth of the 
individual i, with W the total wealth, one has obviously probability that any 
Xi will take a value larger than 1. Therefore, obviously, the distribution of the 
relative wealth Xi cannot be a power law and it has an upper cut-off. Moreover, 
when this cut-off is taken into account the very value of the exponent a (in 
the region of Xi values in which the power law does apply) is affected. 

This discontinuous behavior of the limit A^ ^ oo is important in practical 
applications where the samples are always finite. Indeed, it means that certain 
theoretical results obtained within the assumption of A^ = oo are invalid. The 
problem of finite-size effects in stock market models has been investigated [2], 
especially for the price changes. Here we concentrate mainly on the effect of 
finite-size cut-off on the power-law of the wealth distribution. 

In this paper we highlight the qualitative difference between the behaviors 
of the finite A^ system and the theoretical results for A^ = oo. We consider 
first the example of a pure multiplicative stochastic process that leads for 
A^ = oo to the log-normal distribution. The log-normal distribution implies 
an effective exponent which at asymptotically large values diverges to a — > oo. 
This is compared with the actual behavior for any finite A^ which converges 
towards a = 0. Then, we study the multiplicative stochastic processes with a 
lower bound, which lead theoretically for A^ = cxd to power laws with a > 1 
[3,4]. We contrast it with the finite N systems where values a < 1 can take 
place. 

In particular in financial systems, our results imply that in the absence of a 
social security policy which takes appropriately into account the effects of a 
finite population A^, the Pareto exponent a may fall to values a < 1. This 
in turn means the concentration of an fraction of the total wealth in 

the hands of just a few individuals even if the market size is asymptotically 
large. We show that beyond the moral and social problems which such an 
extreme unequal distribution of wealth implies, a < 1 also means that the 
financial markets themselves become unstable as their fluctuations become 
macroscopic. 
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Let us start with a multiplicative random walk process involving A'^ indepen- 
dent variables (or individual wealths in financial system) Wj, i — 1,...,N, 
which undergo each the updating process: 

Wi(t + l)^Xi(t)wi(t), (2) 



where the random numbers Aj are independent but extracted with the same 
time-fixed and i-independent probability distribution n(A). 

According to the generalized central limit theorem, if the distribution n(A) is 
characterized by the mean 

{lnX)=v, (3) 

and variance 

((lnAf)-(lnA)2 = A (4) 

the asymptotic time dependence of the w probability distribution is: 

P(w, t)dw — , exp ( —(Inw — vtV^ —dw. (5) 

In particular, it was pointed out that this log-normal distribution can be 

regarded as an effective power law (1) with exponent a changing with w [5]. 
For small w, a ~ which corresponds to 1/w distribution, while in the large 
w tail, this distribution (5) behaves as a power law with exponent diverging 
to infinity [6]. 

To see this, let us put the distribution in the form [6]: 

P{w,t)^ 1 exp(aH.O 



where the effective w-dependent power exponent a{w) is defined by 

a{w) = In f . (7) 

^ ^ 2Dt \ex.p{vt)J ^ ' 



Then formally a{w — > oo) = oo, i.e., the tail of the distribution has an effective 
power law with exponent infinity. 

We will show however that this tail is completely irrelevant for any finite 
system at large enough times: there is simply no variable w with such large 
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value within the system. As the time goes on and the distribution spreads 
over larger intervals, the initial number of variables N becomes too sparse 
to sample the far away tail. Thus, in reality, the behavior of the theoretical 
N — oo distribution at asymptotic times is quite different for any system 
containing a finite number of variables. 

An upper bound for the values which the Wi^s can take is given by the total 
wealth in the system, and obviously none of the individuals can have a larger 
value, that is, 

Wi < W, (8) 

with the total wealth W defined as 

TV 

W{t)=Y.Wi{t) = Nw{t). (9) 

i=l 

Thus, from Eqs. (7) and (8) we have 

2Dt \exp{vt) J ^ ' 

In the limit A" = oo, w in Eq. (9) is proportional to the average value 

{w) = JwP{w)dw, (11) 

and for the multiplicative system (2) it is [7] 

(w) = (A)* - exp(ln(A)t). (12) 



However, the result is very different for finite A^. There, the observed mean 
wealth w = J2iWi/N is much less than the A^ = cxo average (w), since the 
average values (and other statistical properties) of the multiplicative system 
are determined by the extreme events and the most extreme events can appear 
only for large enough A^ [8]. It has been shown in Ref. [8] that the mean w 
approaches the value of (w) (Eq. (12)) when the system size A^ > exp(/it) (with 
fi a constant and decided by the details of the multiplicative process) , which at 
asymptotic times is far beyond the size of simulations or real (e.g., financial) 
systems. In fact, for finite-size multiplicative systems the total wealth W{t) — 
J2iWi{t) at large enough times mainly depends on the largest value among Wi 
[7,9], that is, 

W{t) ~ w^^ = ma.x^^^Wi{t). (13) 
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The value of Wmax can be obtained analytically for -C ex.p(Dt), (i.e., for 
finite and asymptotic times) in the case in which the random numbers In A 
are extracted from a Gaussian distribution of zero mean v and finite variance 
D [7,10]: 



~ exp [V2Dt\nN) . (14) 



Substituting Eqs. (13) and (14) into (10) one obtains that 



a<\ . (15) 

\ 2Dt ^ ^ 



Thus, Eq. (15) implies that for finite size N and large enough time t, the 
effective a value is far from being infinite even for the large w tail. In fact, the 
exponent a decreases as time takes larger and larger values. 

The numerical simulations of the system for various times and N values con- 
firm this prediction. One way to analyze the simulation data and to validate 
the power law is to use a Zipf plot [11]. In order to make a Zipf plot one 
re-labels the various values of the individual wealths Wi existing in the sys- 
tem at a given time in descending order as follows. One denotes by w{l) the 
wealth Wmax of the richest individual, w{2) the wealth of the wealthiest among 
the remaining individuals, and so on, until one reaches the poorest individual 
whose wealth is denoted by w{N). With these notations, the power law Eq. 
(1) leads to the Zipf law: 

w{n) ~ n-^/'^. (16) 



Thus, in the log-log plot of w{n) (the Zipf plot) the power law (1) corresponds 
to a straight line with slope —1/a. Visually, steeper slopes represent smaller 
exponents a. E.g., the 1/w distribution (corresponding to a = 0) leads to an 
infinite slope in the Zipf plot, and conversely, a flat and almost horizontal Zipf 
plot corresponds to a — > oo. 

Figs. 1 and 2 exhibit the Zipf plots for the multiphcative system (2) with 
In A distributed by a Gaussian with v = and D = 0.01. As expected from 
Eq. (15), all the slopes for finite A^ are larger than 0, corresponding to power 
laws with effective exponents a much smaller than the theoretical N = oo 
prediction a = oo. In fact, in contrast to the N — oo result, as time increases 
q; — > 0. This is clearly seen in Fig. 1 where the slope —1/a becomes steeper 
with increasing time. In accordance with the same Eq. (15) the Zipf plot (Fig. 
2) gives steeper and steeper slopes {a ^ 0) as N decreases from 5000 to 50. 

The same results can be obtained for different distributions of the random 
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factor In A. Fig. 3 shows the Zipf plots for In A uniformly distributing in a 
range (—0.1,0.1). The above phenomena caused by the upper cut-off (8) are 
still present. 

This upper cut-off effect holds for any finite value of for large enough 
times, and is not limited to the log-normal distribution. A similar argument 
insures that other slow-decaying distributions have sharp truncations and/or 
departures from the N = oo results. In order to uncover these effects one has to 
use appropriate plots for each quantity. For instance, the upper cut-off does not 
affect the straight line shape of the Zipf plot for a system obeying the power- 
law distribution (1). However, this cut-off is made explicit by first normalizing 
each variable (wealth) to the total wealth and averaging over many realizations 
to obtain an estimation of the probability Pareto distribution P{w/W). By 
plotting then log P{w/W) against log{w/W), obviously one cannot expect a 
straight line since at least for w/W > 1 one has log P{w/W) = — oo. In the 
actual plots, indeed a dramatic bent is shown in the graph when the variable 
w approaches the total capital W. 

Thus, we encounter an unexpected subtle effect related with the special prop- 
erties of the power laws. The effect originates in the fact that the density of 
probability P{w) has to be at w = W due to Eq. (8), in spite of the naive 
N = oo expectation [12] that the distribution Eq. (1) will continue to infinity. 

A particular model where these observations are relevant is the multiplicative 
stochastic model with lower bound introduced in Refs. [4,13]. In this model, 
the process (2) is supplemented by a lower bound: 

Wi{t + 1) > cw{t), (17) 

(i.e., whenever according to (2) Wi{t + 1) becomes less than cw{t), the actual 
updated value Wi{t + 1) is taken as cil){t)). Such a lower bound represents in 
the model the economic effects of subsidies, securities, and services. 

For infinite N, this process leads to the precise power law with the exponent 

q; = 1/(1-c). (18) 

In reality, that is, for finite N, one obtains a different result [14]: for c <S 
1/lnN, the exponent is 

Q;^lniV/(lnA^-lnc) < 1, (19) 

while for c 3> 1/lnA'' the relation (18) is recovered. (The large Wi scaling 
properties of this system are representative of a larger class of systems [4,12- 
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17] of the form: 



Wi{t + 1) = \i{t)wi{t) + aiv{t) + b{iv, t)wi{t), 
where b is an arbitrary function [18].) 

Numerical results probing the properties of this multiplicative stochastic model 
with lower bound are shown in Figs. 4, 5 and 6. (We used here the realistic 
asynchronous dynamics, i.e., different wealths are updated in successive ran- 
dom order: at each time t only one wealth is randomly chosen for updating. 
However, similar results are obtained for the synchronous dynamics). In the 
simulations the random factor A is distributed with uniform probability be- 
tween 0.9 and 1.1. 

As explained above, while the Zipf plots of the individual wealth configurations 
do not show any deviation from the precise power- law behavior (sec Fig. 4), 
the Pareto distribution of the normalized w/W, obtained by averaging over 
many runs, shows a sharp bent to log P{w) — — oo for w/W — > 1 (Figs. 5 and 
6). 

In Fig. 5, we have c > 1/lnA^ and then the formula (18) holds especially 
for large N. However, when there is no significant social security pohcy, that 
is, c < 1/lnA^ in the lower bound Eq. (17), one obtains the low values of 
exponent a < 1, as predicted by Eq. (19). This is clearly seen in Fig. 6. 

For the power-law distribution (1), the ratio between the largest wealth and the 
total wealth is Wmax/W ~ A^i/"-i. Thus, when a < 1, the richest person has 
an 0(1) fraction of the total wealth even if the market size is asymptotically 
large. This extreme financial inequality is questionable from moral point of 
view and unacceptable from social point of view [15,16]. 

However, here we make the observation that such an inequality is equally 
dangerous even from the narrower point of view of financial stability. To un- 
derstand the implications of the a < 1 on the market fiuctuations, recall [4,17] 
that the market index fluctuations are expected to also obey (modulo the cut- 
off tails) a power-law behavior with an exponent close to or larger than that 
of the individual wealth Pareto distribution. Telegraphically, the argument in 
Ref. [4] consisted of the following steps: 

• The market index is identified with the total wealth W in the market, i.e. 
with the sum of the individual invested wealths wf, 

• The market fluctuations over a certain time period are therefore the sum 
of the variations of the individual wealths Wi, with i randomly selected for 
updating during that period; 

• The individual variations of the w^'s are, according to Eq. (2), stochastically 
proportional to the WiS themselves, which are distributed by a power law 
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with exponent a; 

• Thus, the market fluctuations are the sums of steps with sizes distributed 
by a power law with exponent a. 

Therefore, a wealth distribution with a < 1 implies market fluctuations dis- 
tributed by a power law with a little larger but still small effective exponent 
(modulo the truncation effects, the details of which are outside the scope of the 
present paper and were studied in Ref. [17]. They fit well the experimentally 
observed cross-over to larger power-law exponents and to exponential-type 
behavior [19,20] for extremely large fluctuations.) In turn this small exponent 
implies the emergence with flnite probability of macroscopically large market 
fluctuations, which means that the flnancial markets in a society with a Parcto 
wealth distribution of a < 1 are open to a signiflcant risk of macroscopic mar- 
ket crashes. The conclusion is that a responsible social pohey which insures 
the smallest wealth Wmin 3> / In A?" is not just a humane duty but also a vital 
interest of the capital markets. 

This finite-size cut-off effect also occurs in systems other than the multiphca- 
tive ones, in particular for the critical phenomena. For the two- or three- 
dimensional percolation system, which has been well studied and understood 
[21], the number Us of clusters containing s sites each decays for large s as a 
power law 

Us ~ l/s^ (20) 

in an infinite lattice and at the percolation threshold Pc- The Fisher exponent 
r is 2.05 for two dimensions (d = 2), while for three dimensions id = 3) it is 
2.19. In a finite lattice, the cut-off behavior similar to Eq. (8) exists, that is, 
no cluster can be larger than the whole lattice and then the cluster radius 

Rs < L, (21) 

in a (i- dimensional lattice of linear size L, where at Pc this radius behaves as 
Rs ~ s^/^, with D = 1.9 for = 2 or 2.5 for d = 3. The size distributions for 
these largest clusters have been investigated [22]. 

Thus, things are more complicated for the cut-off effect of the finite systems. 
For the critical phenomena like the percolation, little is known about the 
extremely rare clusters at Pc with size s between and L'^ which are outside 
the scaling region. As in the case of the multiplicative processes such as Eqs. 
(2)-(17), more work is needed in order to investigate e.g. the margin of the 
wealth value obeying the power-law scaling as well as its dependence on the 
system size. 

Returning to the financial applications, our results caution against the use of 
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the N = oo limit in estimating the parameters of the distributions for wealth 
and market fluctuations. In particular, while for N = oo some random multi- 
plicative mean field models seem to insure a > 1 [12], in reality for any finite 
number of individuals N, in the absence of an appropriate social policy which 
insures the smallest wealth Wmin ^ w/ InN one is led to a very unbalanced 
wealth distribution a < 1. For such a wealth distribution, a large fraction of 
the total wealth remains in the hands of just a few richest individuals. This in 
turn is bound to lead not only to social unrest but also to unbearable financial 
markets fluctuations. 
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Fig. 2. The Zipf plots obtained from the simulations of multiplicative process (2). 
In A is extracted from a Gaussian distribution with v = and D = 0.01. The time 
at which the measurement is performed ist = 10^, and the system sizes are: N = 50 
(+), 500 (*), and 5000 (x). 
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Fig. 3. The numerical Zipf plots of the multiplicative process (2). In A is extracted 
from a uniform probability distribution between the values —0.1 and 0.1. The mea- 
surements are performed at time t = 1.1 x 10^. The system sizes are iV = 50 (+), 
500 (*), and 5000 (x). Inset: N = 500 and different times: t = 10^ (+) and 1.1 x 10^ 

(*)■ 
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Fig. 4. The Zipf plots obtained from the numerical simulations of model (2) with 
the lower bound (17) and asynchronous dynamics, for late time t = 2x 10^ iV and 
different values of c and size N . The random factor A is extracted from an uniform 
probability distribution between 0.9 and 1.1. 
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Fig. 5. Log-log plot of the probability distribution for normalized wealth w/W oi 
the model with lower bound: Eqs. (2) and (17), with simulation parameters: c = 0.3, 
= 100, 1000, and 10000. A is uniformly distributed in the range (0.9,1.1). The 
straight line has the slope 1 + a, with exponent a determined by Eq. (18). 
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